The Effect of a Finite- Width Decision 

■a 

Threshold on Binary Differentially 
Coherent PSK Systems 

By W. M. HUBBARD 

(Manuscript received November 8, 1965) 

A binary PCM regenerator must decide which of two signal states was 
transmitted. This decision is normally based on whether some particular 
voltage is above or below a certain threshold at some particular time. In a 
real regenerator, the voltage in question must differ from this threshold value 
by some finite amount in order to cause the device to respond properly. The 
purpose of this calculation is to examine the effect of this "dead zone" on 
the probability of error in a differentially coherent PSK system both in the 
case where the signal is limited in amplitude (after the noise is added) and 
when it is not. 

It is found that when the expected value of signal power, S, exceeds the 
threshold value, T, by more than 6 db and 10 db in the limited and unlimited 
cases, respectively, the effect on error probability is less than the effect of a 
1-db degradation in signal-to-noise ratio. However, for smaller values of 
signal-to-threshold (S/T) the degradation becomes large. Numerical values 
for error probability as a function of signal-to-noise ratio are presented 
for S/T = oo , 12,9,6, 4, and 3 db for both cases. 

I. INTRODUCTION 

A binary differentially coherent phase shift keyed (DCPSK) system 
is one in which the two (binary) states — "mark" and "space" — are 
transmitted as phase changes between adjacent time slots. In such a sys- 
tem, optimum results are obtained when the expected values of the two 
possible signal states in a given time slot differ in phase by t radians. 
For clarity it will be assumed in the following discussion that the expected 
value of phase change is or t, respectively, according to whether a 
space or a mark was transmitted, even though phase changes of + 
and + t (0 an arbitrary but fixed angle) are equally suitable. 
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Such a system usually consists of a limiter, followed by a storage 
mechanism and a product demodulator, as shown in Fig. 1. If the signals 
in adjacent time slots have amplitudes A and B, respectively, and a 
phase difference ip, the output of the product demodulator is proportional 
to 



v = AB cos \f/. 



(1) 



If the limiter is ideal, A = B = 1. The ideal regenerator samples the 
sign of v and regenerates a mark if v < and a space if v > 0. 

The probability that, due to additive Gaussian noise, v > if a mark 
was sent or equivalently the probability that u < if a space was sent 
has been calculated by many authors. 1 " 4 This probability is given by 



II = fexp (-S/N) 



(2) 



where S/N is the signal-to-noise ratio. 

In a practical system, the regenerator requires some finite value of 
v in order to regenerate a mark or a space reliably. While it is true that 
the effects of this finite-width decision level can be made arbitrarily small 
by providing sufficient gain ahead of the regenerator and sufficient dy- 
namic range for the regenerator, it is often difficult — especially with 
some solid-state devices — to achieve either this gain or the necessary 
dynamic range. For this reason, it may be desirable to operate a regen- 
erator under conditions where the decision threshold is important. In 
this paper we consider the following model : 

(i) An error is made if | v \ > e and, due to additive Gaussian noise, 

v has the wrong sign. 
(ii) No error is made if | v \ > e and v has the correct sign. 
(in) There is a 0.5 probability of error if | v \ < e, regardless of the 
sign of v. 
The probability of error in such a system is clearly dependent on both 
the signal-to-noise ratio S/N and the signal-to-threshold ratio S/T. 
The latter is defined as the ratio of the expected value of signal power 
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Fig. 1 — DCPSK receiver. 
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to the minimum value which causes the output of the product demodu- 
lator, v, to exceed v in magnitude. 

II. SYSTEMS WITHOUT A LIMITER 

If the regenerator were ideal (e = 0), the performance would be inde- 
pendent of whether or not the limiter in Fig. 1 is included, since the sign 
of v is independent of the product AB (both A and B are positive num- 
bers). However, for the case where e > 0, the limiter plays an important 
role under certain conditions. In this section, we consider the case where 
the limiter is omitted. The received signal can then be thought of as the 
vector sum of a unit vector representing the transmitted signal and a 
noise vector with Gaussian-distributed x and y components. The angle 
between the two vectors representing received signals one time slot apart 
together with their magnitudes .4 and B determine v. This is illustrated 
in Fig. 2 for a "mark". From Fig. 2, it is apparent that 

v = AB cos df = AB - A X B X + AyB y 

(3) 

= (A +a x )(B + b x ) +aA 

where A„ = 4-1 or — 1 and B„ — +1 or — 1 depending on the message. 
Cos yp has the wrong sign if the sign of v differs from the sign of A B . 

There will be two probability density functions p~(v) and p+(v), 
respectively, for the cases A B„ = — 1 and A Mo = +1. From symmetry, 
it is apparent that the error probability is the same in these two cases 
The probability of error can thus be determined from either 

II - \ f_ P-(")dv + f p-(v)dv (4) 

or 

IJ = 5 f P+(v)dv + [ p + (v)dv. (5) 

Z J—C J— OB 

Bennett and Salz derive the density functions p±{v) for the quadratic 
form of (3 ) where a x , a u ,b x , and b y are independent Gaussian variables 
of zero mean and variance a 1 . They are given (in our notation) by 

exp[-(.i- 2 + ?/ 2 )/2«r 2 ] 



) 2 + y 2 

(6) 



i r r ex P [-(■■- 
p±(|,) = (2^v .LJ-. VoTlr: 

/ [i»=F (1 +.r)] 2 \ . . 
a 2 is the mean noise power. 
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A ,Bo= UNCORRUPTED SIGNALS 

a,b = noise on A,B respectively 
A,B = received signals 



Fig. 2 — Phasor diagram of a pair of received signals. 
In Appendix A it is shown that substituting (6) into (4) or (5) gives 



+ erfc -je 



— e + ^ 



V2 a y/x* + y 2 



+ i/ 2 
dx dy. 



(7) 



This integral must be computed by numerical methods. The results 
of these computations are shown in Fig. 3. 

III. SYSTEMS WITH AN IDEAL LIMITER 

If an ideal limiter* is included in the system as shown in Fig. 1, (3) 
reduces to 

v = cos ip. (8) 

In this case, it is simpler to perform the calculation of error probability 
from a consideration of the probability density of \f/ itself. From Fig. 4 
and the criteria for error listed in Section I, one obtains 



j /•t+cob -1 (-c) /«coa l e 

TT = i / pO)# + / p(*)# 

2 •'T+coa-ie J x+ c o8-i(-e) 



COS -1 (— C) 



+ir. 

^ •'cos 'e 



(9) 



p(4f)d$ 



* By an ideal limiter is meant a device which removes all amplitude variation 
from the signal without affecting the phase. That is, it transforms the signal 
A(t)e>* u) into the signal A ei* U) . 
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SIGNAL-TO-NOISE RATIO IN DECIBELS 



Fig. 3 — Probability of error in an unlimited DCPSK system. 
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FOR AqAND Bo as shown 
IN FIG. 2, AN ERROR IS MADE 
IF FALLS IN THE DOUBLY 
CROSS-HATCHED REGION AND 
THERE IS A 50% CHANCE OF 
ERROR IF >fi FALLS IN THE 
SINGLY CROSS-HATCHED 
REGION 



Fig. 4 — Geometric interpretation of the limits of integration over pit). 

where p(ip) is the probability density function of ^. The function p(\f/) 
is derived in Appendix B. It is given by: 

p W = -iexp(-y + iexp(^) 

— - — - / cos a cos (a 4- \f/) 

47TO" 2 ■'-(t/2) 

( sin 2 a + sin 2 (a + i/Q \ 
2<r* ) 



•exp 

1 f w " 

-\- - — - / cos a cos (a + ^0 

47TO -2 «'-(t/2) 



da 



r/2 



(10) 



•exp 



( sin 2 a + sin 2 (a + ^) \ 
2^~ ~) 



erf 



cos a 
y/2<r 



cos (a + ft) . 



The probability of error is derived in Appendix C by substituting (10) 
into (9). The result is 



n=-:i 



T/2 



exp 



(1 - e 2 )(V2 
1 - e sin 20 



8 J 



(ID 
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Again the integral must be computed numerically. The results are shown 
in Fig. 5. 



IV. CONCLUSIONS 

From the form of (2) one sees that for the ideal case (no dead zone) 
the probability of error, n(S/X), will plot as a straight line if the scale 
of the abscissa of the graph paper is proportional to the logarithm of the 
logarithm of the scale of the ordinate.* A grid of this type is used in 
Figs. 3 and 5. From these figures one sees that the graph of II vs S/N 
is very nearly linear and very nearly parallel to the graph of the ideal 
case even when the dead zone is significant. 

Since the family of curves of n (S/N) for various values of S/T is a 
set of approximately parallel lines, it is convenient to define a "threshold 
effect noise figure," N T , as the change of S/N which would give a deg- 
radation in error-probability equivalent to that due to the dead zone. 
Since these curves are not exactly parallel lines, N T will be a function 
of the value of II at which it is determined, but this effect is quite small 
over the range of values plotted in Figs. 3 and 5. 

The values of N T as a function of S/T (evaluated at II = 10~ 9 ) are 
plotted in Fig. 6 for both the limited and the unlimited cases. From 
Fig. 6 it is observed that N T ^ 1 db for S/T ^ 9.4 db and 5.8 db in the 
unlimited and limited cases, respectively. However, as S/T becomes 
smaller the effect of N T on error-probability becomes quite important. 

The difference in the value of N T for the limited case and the un- 
limited case gives a measure of the improvement gained from the ideal 
limiting process. This improvement is very small for large S/T; however, 
for S/T = 8 db it already amounts to about 1-db improvement in N T , 
for S/T = 6 db the improvement is about 2 db, and for S/T = 3 db this 
improvement exceeds 4 db. 
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APPENDIX A 

Derivation of the Error-Rate Equation for Unlimited Systems 

Substituting (6) into (4) or (5) and making the change of variable 

* It may be worth mentioning that for 1 db ^ S/N ^ 10 db one can calculate 
the error rate for the ideal case directly on a slide rule by reading 2U on the 
LL3 scale directly opposite S/N (in bels) on the L scale. 
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Fig. 5 — Probability of error in a limited DCPSK system. 
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Fig. 6 — Threshold effpct noise figure vs signal-to-threshold ratio, 
.c —>.c — l gives 
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/ (*-i) 2 + A 
i i r r r exp y 2^ ~) 

3 (27r) 3/ V 2 LJ-J-00 V.r 2 + 1/ 2 



•exp 



+ 



V 2<r 2 [; 



Tri 



•T 2 -f- y 2 } 



V.r 2 + y' 
dx dy dv 



1 



r r r p y 2a 2 ) 

i /_« J_« -vA 2 + w 2 



(12) 



^Tr) 3 ' 2 (T» Je J-oe J-oc 



e **\- 2a'[x> + y*\ ) dxdydV - 
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Making the substitution 



v -\- X 
U = y/2 <r Vx 2 ~+~y 2 



and changing the order of integration gives 

e -|- x 



tt i r r ( u - 1) 2 + A /i rf e + 



z 2 -f w 2 



-5 €rf Va-li + i» + erf Vf7v#TPi d ' 



(13) 



where erf (z) and erfc (z) have the usual definitions 

erf(z) - -7= f exp (-* 2 ) <ft, erfc (a) = 1 - erf (z). (14) 
Vtt •'o 

Equation (13) is readily simplified by means of (14) to give (7). 



APPENDIX B 



(15) 



Derivation of the Probability Density p(i//) 

Bennett 6 and Davenport and Root 6 give the probability density of a 
(see Fig. 2) as 

1 / 1 \ . 1 1 / sin 2 a\ 

p<«) = - exp {- ^j + ^ ^— cos a exp ^- —) 

[cos a "I 

A similar argument gives 

From these results, one can evaluate p(^). 

p(tf) = £ / 8(f + a - P)p(a)p'(0)dadp 

(17) 

= / p(a)p(a + ^)da. 



(16) 
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Cahn 7 gives numerical values for p(\p) for a wide range of values of S/N. 
However, an expression for p(\f/) suitable for calculation of error-proba- 
bilities is obtained as follows. Substituting (15) and (16) into (17) and 
performing the indicated multiplication gives 

p( * } = h L [h exp (" ?) + s v%~* exp ( _ 9?) cos " 

( sin" a\ , con a 1 1 / 1 

• exp { ~ -%r) erf vS + 2 vl^ exp V " ^ 2 

/ sin 2 (a + jgj \ cosJo+jA) 
■CO S (a + *)exp^-- -^ J erf ^ a 

+ -j— i COS a COS (a -f ^) (18) 

4a- 



(sin 2 a + sin 2 (a + ^)\ 
" 2? ) 



cos a; cos (a + \J/) 1 
•erf /— erf 7^ — 71 cos a cos (a 4- i/0 



/ sin 2 a 4- sin 2 (a 4- tfO\\ , 

exp ( 2? ;/*• 



where the other terms vanish because they are odd with respect to 

a — * a 4~ 7T. 

Consider 



/ = 



f ( i\ ( si " 2 (g + *A r cos(« + ^) 
J_ cos (a 4- ^) exp I ^2 ) erf /^ da. (19) 

We observe that the integrand is periodic in period it and the result 
is therefore independent of \f/. Thus, we can write 

f T/2 ( sin 2 a\ cos a 

1 = 2 cos a exp I — , 2 erf /- da 

J-(r/ 2 ) \ 2o- / V2 a 

which is easily integrated by means of the change of variable x = sin a 
and the definition of the error function. The result is 



/ = 2V 2. ff [l-exp(-±)]. 



(20) 



Substituting (20) into (18) and observing that all of the terms of the 
integrand of (18) are even with respect to a — > a -\- t, gives (10). 
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APPENDIX C 

Derivation of the Error-Rate Equation for a Limited System 

Designate the four terms on the right of (10) pi , p 2 , Pa, and p 4 , 
respectively. We observe that 

p<(+) = P4& + *) for i = 1,2,4 

PsGO = -Pstt + *). 

From this we have 



/cos 'e 
(pi(*) +pt(f) +pM)\ty 
- + coa- lE (21) 

/•COS X E 
+ / P 8 (^)#. 

*'»-+oo8 -1 (— e) 

The integration of p x and p 2 is trivial. p 4 can be integrated by changing 
the order of integration and applying (20) twice. Performing the first 
integral in (21) and making a change of variable in the second integral 
gives 

TT = - — — / / ex P ( - o 2 ) dxd y> (22) 

where 

<p = sin -1 e. 

Inspection of the limits of integration in (22) reveals that the integra- 
tion is to be performed over an ellipse centered at the origin. Since the 
integrand is spherically symmetric we can integrate over any quarter of 
the ellipse — say the first quadrant. Then 

COB If 

n-J-»rJP ! ^ s? - p (-5?)'** 

where we have made the change of variable x = r cos 6, y = r sin 0. 
This can be easily integrated over r to give (11). 
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